A sufficient criterion for integrability of stochastic many-body 
dynamics and quantum spin chains 



V. POPKOV 1 ' 4 , M.E. FOULADVAND 2 ' 3 AND G.M. SCHUTZ 1 



1 Institut fur Festkorperforschung, Forschungszentrum Julich, 52425 Julich, Germany 
2 Department of physics, Zanjan university, P.O. box 313, Zanjan, Iran 
Institute for Studies in Theoretical Physics & Mathematics (IPM), P.O.Box 19395-5531, 

Tehran, Iran 

4 Institute for Low Temperature Physics, 310164 Kharkov, Ukraine 



We propose a dynamical matrix product ansatz describing the stochastic 
dynamics of two species of particles with excluded-volume interaction 
and the quantum mechanics of the associated quantum spin chains 
respectively. The time-dependent algebra which is obtained from the 
action of the Markov generator of the exclusion process (or quantum 
Hamiltonian of the spin chain respectively) is given in terms of a set 
of quadratic relations. Analyzing the permutation consistency of the 
induced cubic relations we obtain sufficient conditions on the hopping 
rates (i.e., the quantum mechanical interaction constants) which allow 
us to identify integrable models. From the dynamical algebra we 
construct the quadratic algebra of Zamolodchikov type, associativity 
of which is a Yang Baxter equation. We also obtain directly from the 
dynamical matrix product ansatz the Bethe ansatz equations for the 
spectra of these models. 
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1 Introduction and summary of results 



The notion of integrability in classical stochastic interacting particle systems derives from the 
mapping to quantum spin systems being associated with a solution of the Yang-Baxter equation 



(YBE) [[y, |34[. Classical particle occupation numbers are interpreted as quantum mechanical 
spin degrees of freedom. The generator of the infinitesimal Markovian time evolution thus 
becomes the quantum Hamiltonian of related quantum spin chain [J In its essence integra- 
bility refers to the existence of an infinite set of conservation laws which commute with the 
quantum Hamiltonian and which therefore govern the time evolution of the stochastic process. 
For integrable processes the Bethe ansatz and related methods complement probabilistic ap- 
proaches |J and allow for the derivation of critical exponents, density profiles and correlations, 
shock dynamics and other quantities of interest in the investigation of driven diffusive systems, 
reaction-diffusion processes and other systems both in and and far from equilibrium ||. Thus 
integrability has emerged as a powerful tool in the study of Markov processes involving one or 
more species of particles which move and interact a one-dimensional lattice. Examples include 
the asymmetric exclusion process, spin-relaxation processes and reaction-diffusion systems. 

A major difficulty in making use of integrability is posed by the problem of determining 
whether a given stochastic dynamics actually does correspond to an integrable system since the 
dynamics only determine the quantum Hamiltonian, but does not directly lead to the Yang- 
Baxter equation. The same problem occurs in the context of quantum spin chains alone, i.e., 
without reference to stochastic dynamics. Given the Hamiltonian densities as specified by the 
local interaction constants one would like to know whether the system is integrable or not. 
Even though in special cases an answer can be given by way of straightforward coordinate 
Bethe ansatz || Baxterization [2^] or by an integrability criterion due to Reshetikhin 0, 
a generic method for the construction of the Yang-Baxter equation directly from the quantum 
Hamiltonian is desirable. It is the purpose of this paper to provide such a method for stochastic 
interacting particle systems and their associated quantum spin chain Hamiltonians. 

To this end we extend the dynamical matrix product ansatz (DMPA) j|, |5|], originally devel- 
oped for the single-species exclusion process and the associated spin-1/2 Heisenberg quantum 
chain, to two-species exclusion process with nearest-neighbor interaction and their associated 
spin 1 quantum chains. This is an approach where the action of the time evolution operator 
on an arbitrary state is rephrased in terms of a dynamical algebra of time-dependent operators 
with quadratic relations (see below). In this way we are able to answer two questions. The first 
is: Since the DMPA can be constructed for any stochastic process with nearest-neighbor inter- 
action, what is special about the resulting algebra in integrable cases? We find that the special 
property is associativity which in its turn requires that cubic relations of the operators satisfy 
the YBE. The second question is: What choice of interaction parameters leads to an integrable 
case? Since our method is constructive and yields the YBE without a prior assumption on 
parameters, we obtain relations between the parameters which provide a sufficient criterion for 
integrability directly from the quantum Hamiltonian (i.e., without explicit reference to the infi- 
nite set of conservation laws and to an underlying transfer- matrix) . In this way we recover not 
only a known family of integrable two-species exclusion processes which, using the Reshetikhin 
criterion, was believed to be the only one ||, but also a new family which had been studied 
earlier [|3T], but had not been known to be integrable. 

This paper addresses rather different communities, viz. physicists and mathematicians 



lr The time evolution operators of nonequilibrium systems correspond to non-Hermitian variants of quantum 
spin systems, but in the cases we have in mind this does not affect the integrability. 
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working on interacting particle systems on the one hand and on integrable models on the other 
hand. For the first group some of the remarks on quantum spin chains and integrable quantum 
field theory which appear in the paper may seem obscure and perplexing. Therefore we have 
tried to give a self-contained presentation of our results such that the paper can be followed 
even if all these allusions are being ignored. For a detailed and more pedagogical review we refer 
to ||. Probabilists may also find the mathematically rigorous application of related quantum 
techniques in |ll], |3(J useful in this respect. We hope that for the second group of readers the 
non-standard (and non-rigorous) application of quantum integrability to classical stochastic 
dynamics pursued here will prove inspiring and not distracting. Complementary to the review 
H we refer to [§, [13], [l^, [L4| for the significance of stochastic interacting particle systems in a 
wider mathematical and physical perspective and to [[HJ specifically for driven lattice gases of 
the type studied here. 

The paper is organized as follows. In Sec. II we first follow standard procedure |1 and 
introduce the quantum spin chain representation of the stochastic many-body dynamics. Then 
(Sec. Ill) we extend the DMPA developed in |5| for the single-species exclusion process to 
the two-species case. The special case corresponding to symmetric hopping has been discussed 



recently in p9| . In Sec. IV we derive the YBE from the dynamical algebra and obtain the 
conditions on the parameters which guarantee that the YBE is satisfied. In Sec. V we derive 
the nested Bethe ansatz equations for the spectrum of the quantum Hamiltonian. The idea of 
proposing matrix states as special eigenstates of periodic quantum systems was first employed 

m mil. 



2 Quantum Hamiltonian formalism 

The method of defining stochastic dynamics in terms of an imaginary-time Schrodinger equation 
with a "quantum Hamiltonian" as generator of the Markovian stochastic dynamics has been 
reviewed in detail in ||. In order to be self-contained we present here a summary of what is 
required below. 

Consider a configuration r\ of particles on a lattice of L sites. Here r\ = {r](l), "^(2), . . . , t]{L)} 
where r] (x) denotes the state of the system at site x in terms of an occupation number for 
each permissible particle species. For definiteness we consider here only two-species exclusion 
processes where each lattice site is either occupied with a particle of species A (local occupation 
number = 1) or species B (rzf = 1) or vacant respectively (n^ = = 0). We may therefore 
define the "spin" variable i](x) = — € {1,0,-1} as a unique variable specifying site x. 
Physically one may interprete r](x) e.g. as the charge of the particle occupying site x. 

In the course of time random events take place which change a configuration rj of the system 
with rate w v ^ v '. We define a stochastic process in terms of a master equation 

j t P( V ; t) = J2 K-n, W; t) - w v ^P( m t)} (i) 

v'ex 

for the probability P(r/; t) of finding the state rj at time t. 

The idea of the "quantum Hamiltonian" formalism is to represent each of the possible 
particle configurations i] by a vector | r] ) which together with the transposed vectors ( rj | form 
a basis of the vector space X = (C 2 )® 1 and its dual respectively with scalar product ( r] \ r]' ) = 
5 v rf . We represent an empty site by the symbol and occupied sites by A, B and choose as 
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local basis vectors 




A state r\ of the entire system is then represented by a tensor state 1 77 ) = | r\\ ) <g) . . . <g) . . . | t)l )■ 
Therefore the probability distribution is given by a state vector 

\P(t)) = J2 P (v,t)\v)- (3) 

In this formalism one rewrites the master equation in the form of a Schrodinger equation in 
imaginary time, 

lp( m t) = -( V \H\P(t)) : (4) 

where the off-diagonal matrix elements of H are the (negative) transition rates wi^^u between 
states and the diagonal entries are the inverse of the exponentially distributed life times of the 
states, i.e., the sum of all outgoig transition rates w v ^ v i from state r\. In quantum mechanical 
interpretation r\{x) may be regarded as the ^-component of the spin of an atom in a spin 1 
state. 

The master equation is linear in time and therefore has a formally simple solution. A state 
at time t = t + r is given in terms of an initial state at time t by 

\P(to + T))=e- Hr \P(to)). (5) 

The expectation value Pk(t) = (s\nf\ P(t)), Z — A,B, for the Z-species density at site x is 
given by the projection operator nf which has value 1 if there is a particle of type Z at site k 
and otherwise. The constant vector ( s | = J2 V&X ( rj | performs the average over all possible 
final states of the stochastic time evolution. The real part of the spectrum of H is the set of all 
inverse relaxation times of the process. A non- vanishing imaginary part signals the presence of 
currents characterizing an non-equilibrium system. An equilibrium system satisfying detailed 
balance has a purely real relaxation spectrum and can always be mapped to a quantum system 
with real symmetric Hamiltonian H . 

Here we consider diffusive systems with only hopping processes between neighboring sites. 
We assign hopping rates as follows 



AO - 


-> OA 


with rate gAo 


OA - 


-> AO 


90A 


BO - 


■+ OB 


9bo 


OB - 


■+ BO 


9ob 


AB - 


-> BA 


9ab 


BA - 


-> AB 


9ba 



(6) 
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The transition matrix for the process defined on two lattice sites takes the form 
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(7) 



In a periodic system with L sites one therefore has 

L 

#=5> ( 8 ) 

i=i 

where hi = lcg)...(g)/iCg)...®l is the hopping matrix (|7|) acting on sites i, i + 1 and 1 is the 3x3 
unit matrix. For systems with open boundaries where particles are exchanged with external 
reservoirs the hopping matrix hi is replaced by suitably chosen boundary matrices &i, b^. 



3 Dynamical Matrix Product Ansatz 

In order to introduce the dynamical matrix product ansatz for a three-state system we define 
matrix valued vectors for a single site 



A) = D A (9) 




X) = X A (10) 



with time-dependent matrices E, D A , D B and the time- dependent auxiliary matrices X°' A ' B . 
By taking a L-fold tensor product one obtains a matrix product state (MPS) 



E \* L - 1 


( EQ 


D A ® 


D A Q 


D B I 





P(t)) = I D I (g) D Q . (11) 



These matrices are chosen such that the MPS satisfies the master equation @ for the 
hopping process defined by (0). Notice that there is the multiplication by some additional 
matrix Q of the L-th term in (|TT|). Choosing Q to be time-independent, Q = 0, and following 
01 this leads to an infinite-dimensional algebra of the matrices introduced above and their time 
derivatives which has quadratic relations given by 




|X) ® | A) - | A) ® |X). (12) 
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Special care has to be given to the action of terms h Ll h L _i of the Hamiltonian (H) on P(t), 
due to matrix Q involved in (|11|). It can be shown, however, |2(J that with Q satisfying 

[Q,E- 1 D z ] = [Q 1 D z E- 1 } = 0; Z = A,B. (13) 

(we assume further invertibility of E), no extra relations arise in addition to (|l|). Note that in 
general we assume that Q commutes neither with D nor with E separately, [Q, E] ^ 0; [Q, D] ^ 
0. 

Eqs. (|i~2"|) contain nine relations, namely: 



- (EE + EE) 
~ [ED A + ED A ) + g 0A ED A - g A0 D A E 



2 V 



n + ED B ) + g 0B ED B - g m D B E 



~ (b A E + D A E) - g 0A ED A + g A0 D A E 



- [ d a d a + d a d a 



2 

D A D B + D A D B ) + g AB D A D B - g BA D B D A 



l - [b B E + D B E} - g 0B ED B + g m D B E 
D B D A + D B D A ) - g AB D A D B - g BA D B D A 



I (£>B D B + Y)B Y)B 



X°E- 


EX 


(14) 


X U D A 


- EX A 


(15) 


X°D B 


— EX B 


(16) 


X A E — 


D A X° 


(17) 


X A D A 


— D A X A 


(18) 


X A D B 


— D A X B 


(19) 


X B E- 


D B X° 


(20) 


X B D A 


— D B X A 


(21) 


X B D B 


-D B X B 


(22) 



If one finds matrices satisfying these relations one can calculate expectation values and 
probabilities P(rj,t) by taking a trace in the space on which the matrices act, i.e. 

\P(t)) = Tr\P(t))/Z L (23) 

Here Zl = Tr C L Q with C = E + D A + D B . This is the sum of all unnormalized configurational 
probabilities and hence yields the correct normalization factor. For a given configuration r\ the 
probability P(j],t) is therefore obtained by taking the trace over a suitably chosen normalized 
product of L matrices D A , D B , E. One represents an occupied (vacant) site by a time- dependent 
matrix D z (E) in a string DDDEDEE ... of L such matrices. For illustration consider the 
probability (n z riy {t)) of finding two particles of species Z, Z' at sites x, y in an otherwise 
empty system. One has 

(n z (t)n z '(t) ) = Tr (E x ~ l D z E y - X ~ l D z ' E L - y Q) / Z L . (24) 

The initial state of the system is encoded in the initial values at t = of the matrices. This 
time-dependent algebra generalizes the stationary two-species case studied in detail in |10| and 
the time-dependent single-species algebra introduced in Q. 

The algebra can be exploited either by studying explicit matrix representations or on a 
purely algebraic level [p3[| . Here we choose the second approach which has two objectives: 
(i) the elimination of the auxiliary operators in such a manner that information necessary 
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to calculate the spectrum of H as well as the configurational probabilities is retained, (ii) the 
elimination of the time-dependence from the algebraic relations. Our analysis consists of several 
steps, each of which leads to a reduction of the full algebra ([14]) - (p2|) to a smaller algebra with 
fewer generators and fewer relations. 

Step 1: 

Choosing 

E = (25) 
one can without loss of generality satisfy ([14]) with the choice 

X° = (26) 

One has the remaining eight relations 

^ED A ^+g 0A ED A -g M D A E 
l -(ED B )+g QB ED B -g m D B E 
^{D A E)-g 0A ED A + g M D A E 
I (D A D A + D A D A ^j 
\ [b A D B + D A D B ) + g AB D A D B - g BA D B D A 
^{b B E)-g m ED B + g m D B E 
\ (b B D A + D B D A ) - g AB D A D B - g BA D B D A 

I (pB D B + D B DB ^ 

A similar constraint was used in Ref. [28] for the study of the single-species case, leaving three 
instead of four relations. In |29| the choice C = is made for the two-species case. In this 
case the sum of the auxiliary matrices is set to zero which has been shown to involve no loss of 
generality |33 . 



—EX A 




(27) 


—EX B 




(28) 


X A E 




(29) 


X A D A — 


D A X A 


(30) 


X A D B — 


D A X B 


(31) 


X B E 




(32) 


X B D A - 


D B X A 


(33) 


X B D B - 


D B X B 


(34) 



Step 2: 

Adapting the strategy of Ref. || we further assume that E is invertible. This allows us to 
express the remaining auxiliary matrices X AjB in terms of the physical matrices E, D A , D B . 
Multiplying eqs. ([??]), (^8|) , (^9|) , (|32f) from the right and left resp. with E~ l we find: 



2X A = g A0 (D A + E- 1 D A E)-g QA (D A + ED A E- 1 ) (35) 
2X B = g B0 (D B + E- 1 D B E)-g 0B (D B + ED B E- 1 ) (36) 

Now there are six relations left. Two involve the time-derivatives of D A and D B respectively: 

t) A = -g A0 (D A -E- 1 D A E)-g 0A (D A -ED A E- 1 ) (37) 
t) B = -g B0 (D B -E- 1 D B E)-g 0B (D B -ED B E- 1 ) (38) 
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Using also this yields four other equations: 



g A0 D A E~ l D A E + g 0A ED A E- 1 D A 
g B0 D B E- 1 D B E + g 0B ED B E- 1 D B 
g 0A ED A E' l D B + g m D A E' x D B E 
g 0B ED B E- 1 D A + g A0 D B E- 1 D A E 



lA\2 



(9ao + 9oa)(D j 
( 9 bo + 9ob)(D b ) 2 
(9ao + 9ob-9ab)D a D b 



9baD b D a 



(9bo + g A - g BA )D B D A + g AB D A D 



(39) 
(40) 
(41) 
(42) 



Thus the originally quadratic problem of the time-evolution with nine relations for three physical 
operators and three auxiliary operators has been converted into the quartic problem fl3"7|) - (|42| ) 
with six relations for the physical operators alone which have to be satisfied at all times. Notice 
that the relations fl37|) , fl38|) are linear in the .D-operators while the relations (|39"D - ([42] ) are 
bilinear. 

Step 3: 

The linear relations fl37|), (|38|) are sufficient to describe the one-particle sector, i.e. the set of 
configurations with only one particle on the lattice. In order to eliminate the time-dependence 
we formally define "Fourier" components 



(43) 



D A (t) = a k E k - 1 D A (t)E- k ; D B (t) = f3 k E k ~ L D A (t)E 



k rpk—1 7-1^4 / 



which have the property 

E- l V A (t)E = ae ip V A (t) 
E~ l V B (t)E = f3e ip V B (t) 

Additionally we require that the matrix Q from (|TTD,(|2"4D obeys 

[Q,Dl(t)] = o 

which contains ( |i~3l) special case and yields 

[Q,vZ(t)] = o 



Conversely one has 



Di' B (t) = J- T e-^ B (t)dp. 



The Fourier ansatz turns the time- dependent relations (|37| ) 
order differential equations 

V^ B (t) = -e A ' B V A > B (t) 

with the "dispersion relations" 



g 0A a : e w + g A0 ae ip - g 0A - g A0 
goBf3~ 1 e~ lp + 9Bope w - g 0B - gBo- 



(44) 

(45) 
(46) 

(47) 

(48) 



(49) 

|) into two ordinary first- 

(50) 

(51) 
(52) 
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In this way one can express the time-dependent matrix T>p B (t) in terms of its initial value as 

Vf B (t) = e-^' Bt V^ B (0). (53) 



In what follows we shall omit the time-argument in the initial matrices Drj(0). 



In terms of the Fourier components (|43| ) the four relations ( |39"D - (^) turn into double- 
integral relations where the time-dependence is shuffled into an exponential. Because of (]45l ) 
the quartic relations turn into quadratic relations. Using ( |4~9| ) with pi for the A species and p 2 
for the B species in the third relation and vice versa in the fourth relation below one gets 



= J d Pl J dp 2 a 12 V^ 2 e-^i + ^ (54) 

= J d Pl J dp 2 b 12 V B V B 2 e' { < +t ^ (55) 

= |d Pi y"dp 2 [c 12 P^ 1 P^-^a-V^P^P p A 1 ]e-(4 + ^) i (56) 

= y"dp i y'dp 2 [eZ 12 X>^ 2 ^^ (57) 
with the functions 

a 12 = a(p u p 2 ) = g 0A a~ 2 e- ipi ~ ip2 + g A0 ~ (9oa + g A Q)a~ l e~ ip2 (58) 

612 = b( P i,p 2 ) = g 0B f3- 2 e-^- ip2 + g B0 - (g 0B + ^/TV* 3 (59) 

C12 = c( PuP2 ) = goAa' 1 ^-^-^ + g B0 - {g 0B + g A0 - g AB )r X ^ lV2 (60) 

d 12 = d(p h p 2 ) = ffoBt*- 1 /?-^-* 1 -* 3 + g A0 - (g 0A + 9bo - g BA )cT x z- ip2 (61) 

The four integral equations (0) - fl5"T|), obtained from (p5|) together with the assumption of 
existence of and with the choice (p6| ), form the basis of the subsequent analysis. As an 
intermediate summary we remark that at this point the only relations (out of originally nine) 
that remain are four bilinear relations (0) - (57|). Correspondingly, all expectation values can 
be calculated from the initial matrices using (^), (pT^), (|5^). E.g. one has (see Q23|),(fP|),(P5|)) 

(nf(t)nf (t)) = Tr {D* E L Q) / Z L = J J d Pl dp 2 e-^ + ^> e~ tp ^™Tr (V z p J) z p ' 2 E L Q) / Z L . 

(62) 

The only unknown quantities are time-independent matrix product elements of the form Tr (Dp (0)Dp 2 (0) . . . £ 
to be discussed below. 



Step 4 : 

Before proceeding further a distinction between two different cases must be made. Since we 
assume that no particle species is completely immobile (i.e. the possibilities g A0 = g 0A = or 
9m = 9o B = are excluded) the functions a± 2 and b± 2 do not vanish identically. Indeed in the 
generic case, in the following referred to as case I, none of the four integrands (|5*^ ) - ( |5"7D vanish 
identically. Only for the special case = g A = 9m = g BA = g AB — g A o — go B corresponding to 
C12 = 9 BA = (or the equivalent case obtained by interchanging A and B species) the integrand 
in ( |56"D (or (0) resp.) is zero. This is case II, to be treated separately. 

Case I: 
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The relations (54) - ( |57|) may be reformulated by splitting the integral into two parts as 



pi 



^Pl,P2 d Pl d P2 = / / •••+/ /••• = J P2<pl + Ip2>pl (63) 

TV J —TV J —TV J —TV J —TV J pi 

By changing the order of integration and interchanging pi <-> p 2 in the last term in (|63|) we 
obtain 

/7T rP2 j*TV ppi 

\ Fp^pAV2&P\ = / ^p 2 ,Pl d Pl d P2 (64) 

-7T J —TV J —TV J —TV 



Using fl64|j6"3|) the relation (|54|) thus becomes 



/•pi 

dpi / dp 2 [a 12 Z>£p£ + a^P^e^^* = 0. (65) 



In order to satisfy this for all times i we require the integrand inside the brackets to vanish. 
This is a sufficient condition for satisfying also the original equations fl5"4|). A similar condition 
is obtained for (|5q) . Relations (|56"D , ([57]) can also be satisfied in this manner, but one has to 
require 

+ = + ( 66 ) 

This implies the constraints 



ugAo = (39B0 



on the hopping rates. 

In what follows we assume 



(67) 
(68) 



and (|6q) to hold. One obtains the following four relations 



P\ P2 

B V B 

Pi Pi 

g BA a-\-^V B V2 V^ 

B t^A „ q— 1„— ipi t->A t\B 

Pi 



A 

I'l " P2 

B 

Pi ~ P2 



(69) 
(70) 
(71) 
(72) 



We stress that in these relations the time-dependence drops out. Therefore these are four static 
relations on the operators £>p (0) which together with the two equations (|3|) and with (p5|) , 
( |26| ) form a sufficient (but not necessary) condition for satisfying the original algebra ( |I2| ) on 
the manifold defined by (|67|) , (|68|) . 

These algebraic relations can be written in a compact form by defining the operator valued 
two-component vector 



and the 4x4 matrix 



a 



£(pi,p 2 ) 



AA 
AA 







v° 




) 








AB 
a AB 

AB 
°BA 


BA 
°AB 

BA 
°BA 









a 







BB 
BB 



(73) 



(74) 
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with 

-a = (75) 

ab = a/3c 2 id 12 - gAB9BAe~ ipi ' tp2 . . 

^ a/5c 12 rf 12 - ^ B ^ BA e- 2 ^ 

^AB = 9 BAP n 1 ('') 

a/3ci 2 di 2 - g A B9BAe 2tp2 

ab e-^c 12 -e-^c 21 

^ba = £/Asa— — ^— (78) 

a(3c 12 d l2 - gAB9BAe 1%P2 

ba af3d 21 c 12 - 9ab9ba^ ip ^ IP2 ^ 

° BA ~ apc 12 d 12 - 9AB9BAZ- 2ip2 

OH = (80) 

With these quantities relations (|69|) - ([F^) read 

P P1 g> P P2 = V(p u p 2 )V m <g> V P1 (81) 

We remark that the Fourier components T>p(0) satisfy the algebra of creation operators in a 
1 + 1-dimensional integrable quantum field theory PTfl . The matrix S = SP with the permu- 
tation operator P acting on the two vector spaces may then be regarded as scattering matrix 
with matrix elements Syyi in row YY' and column ZZ', 

= SyKK (82) 



Applying (81) twice shows that S to satisfy 

E(pi,p 2 )S(p 2j j)i) = 1 (83) 

which is fulfilled without further constraints on the hopping rates. This is the analog of the 
field theoretical unitarity condition. 

Case II: 

Case II corresponds to the manifold 

9oa = gso = 9ba = 0, 9ab = 9ao + 9ob (84) 
which is the totally asymmetric two-species exclusion process 

AO — > OA with rate gAo 

OB -> BO g 0B (85) 

AB -> BA g A o + goB- 



investigated in ||31|| . By exchanging vacancies with B particles this process is equivalent to the 
totally asymmetric process with both species hopping to the right, but 

AO — > OA with rate g^o + goB 

BO -> OB g 0B (86) 

AB -> BA g AQ . 
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On the manifold ([84]) relation ( pop is satisfied identically. Without further constraints on the 
independent hopping rates g Ao , g B and on the free parameters a, (3 one may satisfy (0), ( |55| ) 
by splitting the integrals into two domains as in case I, but by requiring the integrand in (|57| ) 
to vanish without splitting of the integral. As in case I the resulting algebraic relations take the 
form (|82[) with 

SH = _*1 = _£Z£^. (87) 
S B A i = S^ B A = (88) 

SiE = ^d 2l/9AB = "toAo + goBe^ 

agABe IP2 

sEj = fr x **9AB/du = Rn a f f ^ oo) 

oBJ? &M e-*i (g-£»») . 

Notice that the corresponding S-matrix is diagonal. The unitarity condition fl83| ) holds. 
The crossing-symmetry relation is generally not satisfied. 



4 Dynamical algebra and Yang-Baxter equations 

The quadratic relations discussed above are sufficient to describe only the sectors with one or 
two particles respectively. In order to study the n-body problem within this approach one has to 
make sure that the expectation values expressed in terms of the quantities Tr (D pi D p ^ D pz . . . E L )/Z L 
automatically satisfy the original master equation, irrespective of the contours of integration 
over the pseudo momenta p^ This can only be ensured by requiring associativity of the algebra 
defined by fl8ll) which in turn implies conditions on the properties of the E-matrix (|74|) or the 
S'-matrix respectively. 

The matrix E acts like a generalized permutation operator on the tensor product of vec- 
tor spaces defined by (|75|). Hence associativity implies that different orders of permutations 
must lead to the same final result. This in not trivial since the vector components are non- 
commutative objects. Let us define SW(p 1 ,p 2 ) = E(pi,p 2 ) <S> 1 with the 2x2 identity matrix 
1 acting trivially on the third subspace of a tensor vector V pi ® V p2 ® V pz . Analogously we 



define E^ 2 ) as acting trivially on the first subspace. Applying (|8~lD in the order 1 
1 <-> 2 yields 



2, 2 



V Pl ®V P2 ®V p . A = ^ip 1 ,p 2 )E^(p 1 ,p 3 )^(p 2 ,p 3 )V P3 



v P2 ®v pl 



(92) 



V2 



pi ■ 



(93) 



On the other hand, choosing the order of permutations as 2 <-> 3, 1 2, 2 <-> 3 one arrives at 
V Pl ®V P2 ®V p , = E^(p 2 ,p 3 )^ 1 \p 1 ,p 3 )^ 2 \p 1 ,p 2 )V P3 G V 
Associativity therefore implies 

S( 1 )(p 1 ,p 2 )S( 2 )(p 1 ,p3)S( 1 )(p2,P3) = S( 2 )(p 2 ,p 3 )S( 1 )(pi,P3)S( 2 )(pi,P 2 ) 



(94) 

This is set of 64 equations for the hopping rates which must be satisfied for all pi. (Extra 
solutions for special values are discussed below). In terms of the elements of the S'-matrix the 
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relations (|94D read 

Sf/' (pi , P2)Si,% (pi , Pz)S 3 -,?l„ (p2 , p 3 ) = 5j^' fc " (p 2 , p 3 ) ' (p a , Ps)^-,, (pi , p 2 ) . (95) 

These are the Yang-Baxter equations with the usual Einstein convention of summing over 
internal indices. In another compact form they may be written 

S12S13S23 = S23S13S12 (96) 

where SV,- is the S- matrix acting on spaces i,j as a function of the pseudomomenta Pi,Pj. If 
these equations are satisfied no extra constraints arise from consistency relations involving more 
than three operators T>. 

Since in case II the ^-matrix is diagonal the YBE is satisfied automatically, i.e., there is 
no further constraint on the hopping rates which would be required for integrability. In case I 
some more discussion is necessary. The various conditions on the rates arising from the YBE 
can be obtained analytically using the software packages mathematica or maple. Due to the 



fact that there is a charge conservation ( S\- = 0, unless % + j = i' + j' ), only the equations 
(|95 ) with i + j + k = %' + j' + k' are nonzero, which leaves 20 equations. Out of them 6 



equations are satisfied trivially. The remaining 14 equations are pairwise equivalent, leaving 
only 7 independent relations. 

It is convenient to pick one of the those independent equations, derive the arising constraints 
on the rates (if any), use this constraint in the remaining equations and then to iterate until 
all equations are satisfied. The solutions that we have found can be classified according to the 
values of the hopping rates gzz'\ = 0, A, B. The crossing-symmetry relation (written as 

S^p (p, q) = Sp^(—p, —q), A = B,B = Afor the spectral-parameter dependent S matrix ( |102| ) 



) is in general not satisfied. In that respect the algebra (82) is not exactly the Zamolodchikov 



algebra ||27||, but the one with some field-theoretical restrictions relaxed. 



4.1 All hopping rates nonzero 

This group consist of a solution with 

(a) a = /3; g A o = 9bo = 9ab = g; 9oa = goB = gBA = h (97) 
and those obtained by relabeling the particles/holes, A <-> B, B <-> 0: 

(b) a = (3; g A0 = g B o = gBA = g; goA = goB = gAB = h (98) 

(c)a = (3h/g; g Ao = g 0B = g AB = g; g 0A = g B0 = g BA = h (99) 

Note that other reshuffling of labels will not result in new sets of rates, e.g. relabeling A «-» 
in (^) gives again (|97|) . We have checked that all the cases (|97| - |99D lead to the ^-matrix of the 
same type given below. Therefore only the solution ( p7|) will be considered in detail. 

Note that constants a, (3 are defined up to common factor, since one can redefine ae tp - 
Therefore in case ( p7|) one can consider a — (3 — 1 without losing generality. 

The corresponding S'-matrix 



ip 



( 



S(pi,p 2 ) 



^i2-(fe+g)e~'Pi 

K 1 2-(h+g)e~ i P2 











K 12 -(h+g)e-^2 

K 12 -ge- zp i-he-* p 2 

K 12 -(h+g)e-^ 







^i2~ge-'P2-fee~'Pi 

K 12 -(h+g)e-^2 

e~ i Pi —Q- i p 2 

3 K 12 -(h+g)e-*P2 









K 12 -(h+g)e~ 
K 12 -(h+g)e- i P2 



\ 



(100) 



12 



K 12 = he- ip '- tp2 + g 



after the transformation 



-ip 



sinh(A) 
sinh(A + r)) 



(101) 



is parametrized to a difference form S(j>i,p2) = S(Xi — A 2 ). The precise form of S as a function 
of a spectral parameter A p6| S(X) is noteworthy (sh = sinh): 



sh(A + r/) S{X) 



( sh(A - 77) 

e-"sh(A) -e A sh(7?) 



V 








-e" A sh(?7) e' ? sh(A) 









\ 



(102) 



sh(A-r?) / 

The anisotropy parameter plays a special role in the theory of 6-vertex model p6| 



Oil C22 1 C21 Cl2 _ Cl2 C21 
a 1" ^21 ^12 ^21 ^12 1 / \ 

A = 11 - - 21 12 — — = cosh(r7) 

9 /Oil Q22 Q21 C12 V " 

Z V °11 22 21°12 



(103) 



If h = g, the transformation (|101|) is to be substituted with e~ ip = (A + i/2)/(X — i/2), 
to retrieve the well-known rational solution of YBE 5(A) ~ XI + iP, P being a permutation 
operator. 

Note however that for the general case h 7^ g, the solution (|102| ) differs from the usual 
trigonometric one due to the spectral parameter dependence in the adiagonal elements. This 
dependence can be removed however by a similarity transformation, and in addition it plays 
no role in the equations for the spectrum. Note also that the sum of S-matrix elements along 
each column is the same for each column, so that S"-matrix is stochastic. 

Finally, note that the choice of rates (|97|) was listed as an integrable case in II. 



4.2 Some hopping rates zero 



Note that Eq. ( T43|J^ ) imply that both auxiliary constants a, (3 are nonzero. In what follows, we 
require this to hold, a/3 7^ 0. 

1- 9ab = 9ba = 0, gAo = gBo, 9oa = goB, a = p. 



Again as in the case (97), one can choose a — 1. The corresponding S- matrix is proportional 
to a permutation operator, S12 = fuP, defined as P^, = 8iji5i'j. This is a tracer diffusion 
process. One can imagine it as usual exclusion process with particles A and B having identical 
dynamics, but different colors p2||. The proportionality coefficient is 



, _ , f r i _ 9ao - (g 0B + 9Ao)e- ip i + gpBe-™-™ 
9aq ~ (qob + g A o)e- ip2 + g 0B e- ip ^ 

2- gBo = gAo = gBA = 0, g 0A = goB<y/(3, g A B = goB(/3 - a)/(3. 

Note that goB = goA + gAB, and since the constants a, (3 are arbitrary, after relabeling holes 
and particles A «-> we obtain the set of constants considered already in Case II. However 
since it is the different physical system, we shall list it independently here. 

The S-matrix is: 
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S(P1,P2) 



1 



V 












-ae ip 2 







l-ae'P2 



(1- 





1 









(105) 



i-/3e*pi 



The S-matrix above can be viewed as the 5-vertex model. The characteristic parameter A for 
this case is pi ]: 

b(a -e '"■) 



A 



cAA cBB 
AA BB 



cBA cAB 
AB BA 



cBB cBA 
BB BA 



a(b 



-IP2 



(106) 



depends on the value of pi- In the spirit of (|113|) the solution of the problem can be reduced to 
one of finding the eigen- values of the transfer-matrix with site-dependent weights. The essential 
property of the above S{jp\,p2) -matrix is that it cannot be transformed into a form containing 
only the difference of the spectral parameters, unlike the case (|102|) . The corresponding Bethe 
equations are given in the next section. Alternatively, one can proceed by relabeling particles 
and holes A <-> and then using ( |111|) . We assume that a ^ [3 in (|105|) since a = (3 falls into 
special tracer diffusion case (solution 1). 

3- g A = 9ob = gBA = 0, g B0 = g A oa/P, gAB = gAoW ~ a )/P- Tnis model can be 
transformed into the previous one by renaming A <-» B and reversing the direction of particle 
motion g A «-> gAo, etc.. The S'-matrix is of type (|105| ). Consequently the equations for the 
spectrum are analogous to ( 118|) and we shall not separately list them here. 



The cases listed above exhaust the list of the nontrivial solutions of the YBE for two species 
of particles. 



5 Spectral equations 

Because of particle number conservation eigenstates and their eigenvalues can be classified 
according to number N z of particles of each species that move on the ring. We denote by 
iV = N A + N B the total number of particles. The quantity Q = N A — N B shall be referred to 
as charge. 

N = 0: 

This is the empty lattice. Since this is obviously an invariant state under the stochastic dy- 
namics. Hence the single eigenvalue 

e = (107) 

of H in this sector vanishes by construction. 
N =1: 

In order to obtain the relaxation spectrum from the eigenvalues 

4 = gozl-'e- 1 ' + gzol^ p - g oz - g zo (108) 

where 7 = a,j3, Z = A,B depending on the species of a single-particle system on a finite 
lattice, one uses the Fourier ansatz and the solution (^). This requires calculating Tz{p) = 
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Tr (V z (p)E L Q) at time t = 0. Because of the cyclic property of the trace the equation Tz^P) = 
r y L e ipL J-'z(p) must be satisfied. This yields the spectral equation 

j L e ipL = 1 (109) 

and thus fixes the allowed values of the pseudomomenta p = 2nn/L + i In 7. We stress that this 
quantization of the pseudomomentum is not a property of the matrices T>q(p) themselves. It 
appears only as a result of taking the trace in the one-particle sector. For expectation values in 
higher sectors we shall obtain different quantization constraints. We remark that for calculating 
only the spectrum the normalization factor Zl is not required. 

N > 2: 



In the case when there are n particles of both types A and B in the system, the averages 
are written in terms of quantities 



F z ^~ z »{ Pl , ...p N )=Ti {V z \ Vl ) . ..V Z "( PN )E L Q) 
Now we have to distinquish between different integrable models. 



110) 



Case II 

The easiest to handle is the Case II, where S-matrix is diagonal. 

Commuting the T>(pi,) through around the "circle" inside the trace in ( |110| ), and using the 



algebra (|82|) , one gets, using the trace property: 



TV 

n 



s z z k £(PkiPiy pkL 



1, k 



1,- 



. iV 



111^ 



with S'-matrix elements written in ([57|-[n]). a, (3 can be set both to 1 by simple rescaling of 



quasiimpulses of A and B particles. It can be shown [[[7]] that applying the coordinate Bethe 
Ansatz directly to the process (E^) yields the same result (|111|). 



Case I, all nonzero hopping rates 

Commuting the T>(pi) through around the "circle" inside the trace in ( |110|) , and using the 
algebra (|82|), one obtains: 

N-l 

s^(pi,p2) n ^^(pi.p^o 1 ^^'^)^^"^ = f ziZ2 - zn (112) 

i=3 

(summation over repeated indices is implied), that can be shortly rewritten in matrix form as 

- Tr (L 01 (p 1 ,p 1 )L 02 {pi,P2) ■ ■ ■ L m (p 1 ,p N ))J r N = z~ lVkL T N (113) 

using the property S^^p^p) = —5 a p'bpa' of the S'-matrix ( |100| ). The matrix L k{pi,P2) is a 
matrix ( |100j) acting nontrivially in su(2) Q ® su(2)k, and acting as identity matrix in the other 
subspaces from sM^ofT^ ®su(2)j, and Tr denotes trace over the su(2) . 

15 



Subsequent analysis of the above eigenvalue equations can be done in the framework of 
standard coordinate |18| or algebraic nested ]T5[ Bethe Ansatz, leading to following spectral 



equations (we made the transformation ( |101| ) to the difference form): 

sh(X k ) \ L , ,n+\ tt sh(\ k — \ n — v) tt sh(\i^ — A fc — r?) 



1 ^i^ 1 TT 1 fe ~ n ~ n TT y - k - l} (ma\ 



s/i(A fc + 77) / £J s/i(A„ — A fc — 77) a = a ^(Ar - A fe ) 

-«iv TT sh(X^ ~ A n ) = /_i\at s +i TT sh(X ( p - a[/ X) - 77) 
M s/ i (Ai 1) - A n - 77) ' sh(\P - X^ - 77) 



c 



where e 77 = k = 1, 2 ... N, v = 1, 2 . . . N B , N B < N. The above Bethe Ansatz appeared 

without the proof in ||. 

Case I, some zero hopping rates (tracer diffusion) 

Commuting the term depending on (j>k) the "circle" inside the trace in (|110|) , one has using 
the C PL10|J104|) : 



where 

JV 

R k (pi,...p N ) = l[f(Pk,p ] y pkL (117) 

The set of equations contained in (|116|) means that all Rk are strictly the same, R± — R2 — 
. . . R n- Note additionally that in the right-hand side of Eq. ( |116| ), the lower indexes are shifted 
one step to the the right. It gives the guideline for determining the spectral equation. E.g. if 
all Zi = A, shift of the sequence {Z{} leaves it invariant, so we will obtain R^ = 1, for any k. 
If the sequence {Zi} is periodic with period n, recurrent use of ( |116| ) yields i?£ = 1. In the 



general case, when {Zi} is non-periodic, recurrent use of ( 116 ) N times gives: 



Rkipi, ■ ■ -Pn) n = 1, Ri = R 2 = ...R 



Case I, some zero hopping rates (5-vertex model) 

Proceeding analogously to ( [L12| ), and denoting a(p) = 1 — ae ip , b(p) = 1 — f3e w one obtains: 

(« > =m> u- 1 ; S^n a(p , ))a(pt) m 

^ V/3 J ii a(pW)a( P „) ( ' IWW') 1 ' 

where k = l,2...N, v = 1,2... N B , N B < N. 

Alternatively, one can obtain the nested Bethe Anzatz in this and other cases by an addi- 
tional (nested) Fourier transform of either T>a{p) or T>b{p), i.e. 

V B (p) = V A (p) J **»&{q)dq (120) 
16 



We remark that the Zamolodchikov-type algebra not only leads to the spectral equations 
( |114| )-( |TT9| ) but also implies functional relations for the matrix elements ( |110| ). These functional 



relations are satisfied by Bethe wave functions as in the spin 1/2 case |23[ and hence yield 
expressions for the expectations as integrals over appropriately chosen contour. Alternatively 



one could search for representations |35| and directly calculate the matrix elements. This 
procedure requires further investigation. 
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